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.^ ■ Abstract. We consider fast magnetic reconnection that takes place within turbu- 

' lent magnetic flux and show that the process results in diffusion of magnetic fields and 

^D , matter, which we term reconnection diffusion. The process of reconnection diffusion 

is based on the model of 3D reconnection of weakly turbulent magnetic fields and is 
applicable to both fully ionized and partially ionized gas. The rate of reconnection dif- 
fusion does not depend on the level of ionization and therefore the usually employed 
ambipolar diffusion idea gets irrelevant for magnetic field transport in turbulent fluids. 
We claim that the reconnection diffusion process is a manifestation of the violation of 
(-H . flux conservation in highly conducting turbulent fluids. We discuss the consequences of 

f^ ■ reconnection diffusion for star formation and stress. We show that reconnection diffu- 

I I sion on large scales is independent of small scale magnetic field dynamics of magnetic 

O . fields. We conclude that numerical simulations correctly represents the diffusion of ac- 

^ ' tual astrophysical magnetic fields in flows with substantially larger Lundquist numbers 



^ I if these simulated regions regions are turbulent. 



> 

CO' 1. The problem of flux freezing 

m ■ 
m 



Interstellar media are known to be turbulent and magnetized, and both turbulence and 
magnetic field are important for star formation (see Armstrong et al. 1994, Chepurnov 
& Lazarian 2009, Crutcher 2012). The existing star formation paradigm has been de- 
^ . veloped with the concept that flux freezing holds, i.e. that the magnetic field is well 

coupled with ions and electrons in the media (Alfven 1942). Therefore it has been as- 
sumed assumed that the way to move matter across magnetic field lines should involve 
ambipolar diffusion (see Mestel 1965, 1966, Mouschovias et al. 2006). According to 
the theory, the change of the flux to mass ratio happens due to ambipolar diffusion, 
j^ I i.e. to the drift of neutrals which do not feel magnetic fields directly, but only through 

" " ion-neutral collisions. Naturally, in the presence of gravity, neutrals get concentrated 

towards the center of the gravitational potential while magnetic fields resist compres- 
sion and therefore leave the forming protostar. The rate of ambipolar diffusion for a 
cloud in gravitational equilibrium with the magnetic field depends only on the degree 
of ionization of the media. 

Magnetic fields are important at all stages of star formation. In many instances the 
ideas of star formation based exclusively on ambipolar diffusion have been challenged 
by observations (Troland & Heiles 1986, Shu et al. 2006, Crutcher et al. 2009, 2010a, 
see Crutcher 2012 for a review). 

While astrophysical fluids show a wide variety of properties in terms of their colli- 
sionality, degree of ionization, temperature etc., they share a common property, namely, 
most of the fluids are turbulent. The turbulent state of the fluids arises from large 
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Reynolds numbers Re = LV/v, where L is the scale of the flow, V is it velocity and V 
is the viscosity, associated with astrophysical media. Note, that the large magnitude of 
Re is mostly the consequence of the large astrophysical scales L involved as well as the 
fact that (the field-perpendicular) viscosity is constrained by the presence of magnetic 
field. 

The key concept of magnetic flux freezing has been challenged recently. On the 
basis of the model of fast magnetic reconnection in Lazarian & Vishniac (1999, hence- 
forth LV99), Lazarian (2005) claimed that the removal of magnetic fields from turbulent 
plasma can happen during star formation due to magnetic reconnection rather than slow 
ambipolar drift (see also Lazarian & Vishniac 2009). The process that was later termed 
reconnection diffusion does not depend on the degree of ionization, but rather on the 
properties of turbulence. The numerical confirmation of the idea was presented for 
molecular clouds and circumstellar accretion disks in Santos-Lima et al. (2010, 2012). 
The numerical testing of the LV99 reconnection model that underpins the concept of 
reconnection diffusion was successfully tested in Kowal (2009, 2012). Compared to 
this testing of reconnection diffusion, the testing in Kowal et al. (2009, 2012) was per- 
formed with much higher numerical resolution of the reconnection region and much 
better control of turbulence and other input parameters. In addition, more recent formal 
mathematical studies aimed at understanding of magnetic field dynamics in turbulent 
fluids supported the LV99 predictions (see Eyink, Lazarian & Vishniac 2011). 

We note that the interstellar medium is collisional (Yamada et al. 2006), and there- 
fore ideas of coUisionless reconnection (see Shay & Drake 1998, Shay et al. 1998, Bhat- 
tacharjee et al. 2005, Cassak et al. 2006) are not applicable. In addition, a theory of 
magnetic reconnection is necessary to understand whether reconnection is represented 
correctly in numerical simulations. One should keep in mind that reconnection is fast 
in computer simulations due to high numerical diffusivity. Therefore, if there are situa- 
tions where magnetic fields reconnect slowly, numerical simulations do not adequately 
reproduce the astrophysical reality. This means that if coUisionless reconnection is in- 
deed the only way to make reconnection fast, then the numerical simulations of many 
astrophysical processes, including those in interstellar media, which is collisional at the 
relevant scales, are in erroo 

To understand the difference between reconnection in astrophysical situations and 
in numerical simulations, one should recall that the dimensionless combination that 
controls the resistive reconnection rate is the Lundquist numbeo defined as 5 = L^Va/^, 
where L, is the length of the reconnection layer, Va is the Alfven velocity, and X = 
ric^/471 is Ohmic diffusivity. Because of the large astrophysical length-scales L^ in- 
volved, the Lundquist numbers are huge in astrophysics, e.g. for the ISM they are 
about 10^^, while present-day MHD simulations correspond to 5 < 10^. As the numer- 
ical efforts scale as L^, where L^ is the size of the box, it is feasible neither at present 
nor in the foreseeable future to have simulations with realistically Lundquist numbers. 
As the brute approach fails, it is important to understand if and when numerical cor- 



' Note, that it would be wrong to conclude that reconnection must always be fast on the empirical grounds, 
as solar flares require periods of flux accumulation time, which correspond to slow reconnection. 

The magnetic Reynolds number, which is the ratio of the magnetic field decay time to the eddy turnover 
time, is defined using the injection velocity v/ as a characteristic speed instead of the Alfven speed Va, 
which is taken in the Lundquist number. 
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Figure 1. Upper panel. Sweet-Parker reconnection. A is limited by resistivity and 
is small. Lower panel: reconnection according to the LV99 model. A is determined 
by turbulent field wandering and can be large. From Lazarian, Vishniac & Cho 
(2004). 



rectly represent astrophysical reality in spite of the tremendous difference in Lundquist 
numbers. 



2. 3D reconnection in weakly turbulent fluid 

If magnetic fields cannot reconnect fast, the intersection magnetic filed lines should cre- 
ate 3D elastic structure making magnetized fluid Jello-like. This would automatically 
mean that numerical simulations of magnetized fluids cannot represent astrophysical 
phenomena which Lundquist numbers are much larger than those of the simulations. 
To understand whether fluid-like behavior is possible we should consider the process 
of magnetic reconnection. 

The corresponding study of 3D magnetic reconnection in the presence of turbu- 
lence was performed in LV99. According to the latter study, magnetic reconnection of 
a turbulent field is fast, i.e. independent of resistivity and it happens at a rate 

Vrec,LV99-VAiL,/Ly^^Ml for L, < L, (1) 

where L^ is the horizontal extent of the reconnection zone (see Figure [T]), L is the in- 
jection scale of turbulence, Va is the Alfven velocity and Ma = Vl/Va is the magnetic 
Mach number. In Eq. ([Hi it is assumed that turbulence is subAlfvenic. For L^ > L, LV99 
got the expression that differs from Eq. ([T) by a change of the power 1/2 to —1/2. 

The LV99 model of magnetic reconnection generalized the classical Sweet-Parker 
model of reconnection (see upper plot in Figure [1} for the case of turbulent fields. The 
enhancement of the reconnection rate compared to the Sweet- Parker model is achieved 
by increasing the outflow region A by accounting for magnetic field wandering. Indeed, 
in the Sweet-Parker model the reconnection speed is limited by the outflow of the matter 
from the reconnection layer and A is determined by Ohmic diffusion. Therefore, 

A 

Vrec = VA— (2) 
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is much less than Va since A <C L^. At the same time A becomes comparable with L for 
V --^Va provided that A is determined by the turbulent field wandering. 

3. LV99 model and Richardson diffusion 

At the scales less that the scale of injection of the strong MHD turbulence, i.e. for 
I < Itrans (LV99, Lazarian 2006): 

ltrans-L{VL/VAf=LMl. (3) 

the magnetic field lines exhibit accelerated diffusion. Indeed, the separation between 
two particles dl{t)/dt ~ v(Z) for Kolmogorov turbulence is ^ (Xtl^l^, where a is pro- 
portional to a cube-root of the energy cascading rate, i.e. a^ ^ V^/^ for turbulence 
injected with superAlvenic velocity Vl at the scale L. The solution of this equation is 

l{t) = [ll'^ + at{t-to)fl^, (4) 

which provides Richardson diffusion or P '-^ t^. The accelerating character of the pro- 
cess is easy to understand physically. Indeed, the larger the separation between the 
particles, the faster the eddies that carry the particles apart. 

Eyink, Lazarian & Vishniac (2011) showed that the presence of Richardson diffu- 
sion induces fast magnetic reconnection and re-derived expressions for the reconnection 
rate from the Richarson diffusion concept. This provides a way to understand the ap- 
plicability of LV99 approach to more involved cases, e.g. for the case of fluid with 
viscosity much larger than magnetic diffusivity, i.e. the high Prandtl number fluids. A 
partially ionized gas is an example of such a fluids, which is essential to understand in 
terms of star formation. 

In high Prandtl number media the GS95-type turbulent motions decay at the scale 
l±,crit, which is much larger than the scale of at which Ohmic dissipation gets important. 
Thus over a range of scales less than l±^crit magnetic fields preserve their identity and 
are being affected by the shear on the scale l±,cnt- This is the regime of turbulence 
described in Cho, Lazarian & Vishniac (2002) and Lazarian, Vishniac & Cho (2004). 
In view of the findings in Eyink et al. (2011) to establish when magnetic reconnection 
is fast and obeys the LV99 predictions one should establish the range of scales at which 
magnetic fields obey Richardson diffusion. It is easy to see that the transition to the 
Richardson diffusion happens when field lines get separated by the perpendicular scale 
of the critically damped eddies l±^crit- The separation in the perpendicular direction 
starts with the scale r,„,Y follows the Lyapunov exponential growth with the distance I 
measured along the magnetic field lines, i.e. r„„iexp(Z//|| ^hy). where l\\,cnt corresponds 
to critically damped eddies with Iperp.cnt- It seems natural to associate r,-,„v with the 
separation of the field lines arising from Ohmic resistivity on the scale of the critically 
damped eddies 

riit=^h,^rit/VA, (5) 

where r\ is the Ohmic resistivity coefficient. 

In this formulation the problem of magnetic line separation is similar to the anisotropic 
analog of the Rechester & Rosenbluth (1978) problem (see Narayan & Medvedev 2003, 
Lazarian 2006) and therefore distance to be covered along magnetic field lines be- 
fore the lines separate by the distance larger than the perpendicular scale of viscously 
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damped eddies is equal to 

Lrr ~ llcrit ln(/x,cri/ /rinit ) (6) 

Taking into account Eq. ^ and that 

^±,crit=^h\,crit/VA, (7) 

where v is the viscosity coefficient. Thus Eq. Q can be rewritten 

LRR^ll„it'^r\Pt (8) 

where Pt = y/r\ is the Prandtl number. . 

If the current sheets are much longer than LrR, then magnetic field lines undergo 
Richardson diffusion and according to Eyink et al. (2011) the reconnection follows the 
laws established in LV99. In other words, on scales significantly larger than the viscous 
damping scale LV99 reconnection is applicable. At the same time on scales less than 
Lrr magnetic reconnection may be sIowq Somewhat more complex arguments were 
employed in Lazarian et al. (2004) to prove that the reconnection is fast in the partially 
ionized gas. For our further discussion it is important that LV99 model is applicable 
both to fully ionized and partially ionized plasmas. 

4. Reconnection dilTusion model 

The rate given by Eq. ([T|l is sufficient to allow the magnetic field to disentangle during 
the turnover of the eddies. Such eddies are a part of the picture of strong Alfvenic 
turbulence (see Goldreich & Sridhar 1995, henceforth GS95)0 As a result, mixing 
eddy-type motions are possible in high Lundquist number fluids. 

The peculiarity of reconnection diffusion is that it requires nearly parallel mag- 
netic field lines to reconnect, while the textbook description of reconnection frequently 
deals with anti-parallel description of magnetic field lines. One should understand that 
the situation shown in Figure [T] is just a cross section of the magnetic fluxes depicting 
the anti-parallel components of magnetic field. Generically, in 3D reconnection con- 
figurations the sheared component of magnetic field is present. This component is also 
frequently referred to as "guide field" in the reconnection literature. 

As we discuss below reconnection diffusion is closely connected with the recon- 
nection between adjacent Alfvenic eddies (see Figure|2]l. As a result, adjacent flux tubes 
exchange their segments with entrained plasmas and flux tubes of different eddies get 
connected. This process involves eddies of all the sizes along the cascade and ensures 
fast diffusion which has similarities with turbulent diffusion in ordinary hydrodynamic 
flows. 



Incidentally, this can explain the foiTnation of density fluctuations on scales of thousands of Astronomical 
Units, that are observed in the ISM. 

^The Alfvenic cascade develops in compressible MHD turbulence independently from the cascade of 
fast and slow modes, and its properties in compressible and incompressible turbulence are very similar 
(Cho & Lazarian 2003, Kowal & Lazarian 2010). For our arguments related to reconnection diffusion the 
compressible motions are of secondary importance. Indeed, the LV99 reconnection model is governed by 
field wandering induced by the Alfvenic cascade and turbulent mixing and also depends on the solenoidal 
component of the fluid. 
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Figure 2. Reconnection diffusion: exchange of flux with entrained matter. lUus- 
tration of the mixing of matter and magnetic fields due to reconnection as two flux 
tubes of different eddies interact. Only one scale of turbulent motions is shown. In 
real turbulent cascade such interactions proceed at every scale of turbulent motions. 



5. Theoretical expectations 

Ions and neutrals closely follows laminar magnetic field lines. Turbulent reconnection 
requires that magnetic flux gets leaky and plasma move perpendicular to the mean mag- 
netic field. In this situation the mass loading on magnetic field lines changes and the 
mass to flux ratio does not stay constant. 

In the absence of gravity reconnection diffusion would tend to distribute magnetic 
flux and volume uniformly in the volume. Therefore any initial correlation of density 
and magnetic fields should dissipate due to reconnection diffusion. 

In the presence of gravity one should expect density to diffuse towards the grav- 
itating center, while leaving magnetic field behind. This would happen even in purely 
ionized plasmas, as reconnection diffusion does not depend on the degree of ionization. 

Plasma motions along magnetic field lines is expected to be diffusive due to mag- 
netic field wandering. Therefore the concepts of collection of matter along magnetic 
field lines must be modified in the presence of turbulence. 



6. Numerical testing of reconnection diffusion 



Numerical testing of reconnection diffusion concept have been performed in a series 
of papers by Santos-Lima et al. (2010, 2012). Both diffusion in the absence and in 
the presence of gravitational potential was tested. The results corresponded well to 
the theoretical expectations. In particular, in Santos-Lima et al. (2012) the study of 
the formation of circumstellar accretion disks has been performed. In the absence of 
turbulence the disks did not form, demonstrating the well-known magnetic breaking 
catastrophe behavior. Indeed, ambient interstellar- magnetic field were getting stronger 
during the phase of accretion and were very efficient in removing angular momentum. 
Reconnection diffusion, however, was removing magnetic fields from the turbulent ac- 
creting matter and was shown to be able to resolve the magnetic breaking paradox. This 
corresponded to the theoretical expectations in Lazarian & Vishniac (2009). 
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Numerical modeling of molecular clouds in Santos-Lima et al. (2010) demon- 
strated the ability of matter to concentrate towards the gravitational potential leaving 
magnetic field behind. The effect was shown to be present for different initial condi- 
tions. 



7. Comparing reconnection diffusion predictions and observations 

A comparison of the reconnection diffusion expectations with observations was per- 
formed in Lazarian, Esquivel & Crutcher (2012). In particular, results of Crutcher 
(2010) on the difference of the magnetization of cores and evelopes was analyzed. The 
traditional ambipolar diffusion requires that the less ionized and therefore more diffu- 
sive core should lose magnetic field faster than the envelope that has larger degree of 
ioniziation. There a larger mass to flux ratio is expected in the core compared to the 
envelope. 

Contrary to the above expectations the observations indicated that the mass to flux 
ratio is larger in the envelope than in the core. This cannot be explained with ambipo- 
lar diffusion. At the same time, reconnection diffusion presents us with a number of 
possibilities. As the reconnection diffusion does not depend on the degree of ionization 
this removes the constraint that the magnetic field diffusion should happen faster out of 
the core. For instance, in a model of a homogeneous molecular cloud boadering with 
the diffuse media with lower magnetization, the central part of the cloud keeps larger 
magnetization than its periphery. In the case of the core and the envelope system, the 
core with the lower level of turbulence and larger magnetic field may be subAlfenic, 
i.e. have the Alfven Mach number Ma < 1 and exhibit the reconnection diffusion that 
is suppressed by the factor M^ < 1 . Due to this factor, reconnection diffusion has an 
interesting property, namely, the larger initial magnetization, the slower the diffusion. 
Thus suggests that in the case of a collapse that compresses magnetic field the regions 
of higher magnetization will tend to lose magnetic flux slower, increasing the contrast 
between regions of low and high magnetization. This is provided that the regions have 
the same size. Otherwise, for the same degree of magnetization smaller regions will 
tend to lose the excess of magnetic field quicker. Indeed, the diffusion coefficient K 
scales as l^'^, while the rate of diffusion scales as k/Z^ and therefore r^ l^^'^. In the 
presence of self-gravity, however, one expect the denser regions to decouple from the 
cascade and the turbulent relation to be modified, making cores substantially less tur- 
bulent than envelopes. 

All in all, while the decrease of the core size tends to make the diffusion out of the 
core faster, other effects, i.e. decrease of the turbulence in the core due to decoupling 
of the motions in the core from the outside cascade and smaller Alfven Mach number 
of the core tend to prolong the time of the diffusion of magnetic flux out from the 
core. Due to the effects discussed above, the inhomogeneity of the envelope may also 
contribute to the higher mass to flux ratio observed in the envelope. 

An additional effect related to reconnection diffusion may be important for en- 
hancing the mass to flux ratio in envelops compared to cores. If the matter is accu- 
mulated along magnetic field lines due to external forcing rather than self-gravity, the 
matter will not follow the hour glass morphology of the field, but due to field wandering 
should preferentially reside in the envelope. 

All these processes should be quantified better, but, it is clear that, unlike the 
process of ambipolar diffusion, reconnection diffusion allow situations that the cores 
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are more magnetized than envelops. Further research of the issue as well as more 
extensive comparison of the observational data with model predictions are required. 

Lazarian, Esquivel & Crutcher (2012) compared the rates of reconnection diffu- 
sion of magnetic fields out of the cloud with the rate of gravitational collapse and calcu- 
lated the critical densities at which clouds should collapse trapping the magnetic flux. 
The results of these calculations are encouraging. The column densities of the clouds 
at which the magnetic trapping happens were found to correspond to observations in 
Crutcher (2012). 

It is encouraging that the reconnection diffusion concept provides new ways of 
approaching challenging problems of star formation in extreme environments. For in- 
stance, galaxies emitting more than 10^^ solar luminosities in the far-infrared are called 
ultra-luminous infrared galaxies or ULIRGs. The physical conditions in such galaxies 
are extreme with a very high density of cosmic rays (see Papadopoulos et al. 201 1). 
Ambipolar diffusion is expected to be suppressed due to cosmic ray ionization. At the 
same time these environments have the highest star formation rates. This is suggestive 
of a process that removes magnetic fields independent of the level of ionization. Re- 
connection diffusion does not depend on ionization and provides a natural explanation. 

It looks that reconnection diffusion can provide solutions to other long standing 
problems. For instance, the observed star formation rate is about the same in galaxies 
with low metallicities as in galaxies with high metallicities (see Elmegreen & Scalo 
2004). This is difficult to understand if magnetic field loss is governed by ambipolar 
diffusion, that is supposed to be much faster in low-Z galaxies with high metallici- 
ties. As a result, ambipolar diffusion theory predicts that the Initial Mass Function 
(IMF) and therefore the star formation rate would be shifted in low-Z galaxies. This 
prediction contradicts to observations. As we discussed earlier, reconnection diffusion 
efficiency does not depend on metallicity, and this provides a solution consistent with 
observations. 

In addition, an important property of reconnection diffusion is that, unlike ambipo- 
lar diffusion, it not only provides the removal of the magnetic field, but also induces its 
turbulent mixing, which tends to make the distribution of magnetic field uniform. Thus 
turbulent mixing mediated by reconnection diffusion helps to keep the diffuse media in 
a magnetized subcritical state. In this situation, the external pressure is important for 
initiating collapse, which well corresponds to the observations of numerous small dark 
clouds not forming stars in the inter-arm regions of galaxies (Elmegreen 2011). 

Finally, the reconnection diffusion concept implies that in realistic turbulent media 
there is no characteristic density for the collapse to be initiated. Therefore any cloud 
with the appropriate virial parameter (see McKee & Zweibel 1992) can form stars. The 
difference between different clouds arises from the density controlling the timescale of 
the collapse and kinetic energy determined either by the level of steering by turbulence 
or thermal motions of particles. Directly, the requirement of clouds to be molecular is 
not present for the reconnection diffusion to induce star formation. The only difference 
between molecular and atomic clouds is that the former have lower temperatures. 



8. Reconnection diffusion and MHD numerical simulations 

Reconnection diffusion is based on the numerically tested LV99 model of fast reconnec- 
tion in turbulent media as well as more recent insights into the violation of the frozen-in 
condition in turbulent magnetized fluids (see Eyink et al. 2011). The direct 3D MHD 
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simulations of the magnetic field diffusion in Santos-Lima et al. (2010, 2012) showed 
consistency with the theoretical expectations, but such simulations on their own cannot 
be used to justify the concept of reconnection diffusion. The interpretation of the results 
requires the proper understanding of scaling of magnetic reconnection with the dimen- 
sionless combination called the Lunquist number S = {LcuKshVA/'f]), where Lcur.sh. is 
the extent of the relevant current sheet, r\ is Ohmic diffusivity. The Lundquist numbers 
in molecular clouds and in the corresponding simulations differ by a factor larger than 
10^. In this situation one can establish the correspondence between the numerical sim- 
ulations and astrophysical reality only if the reconnection does not depend on S. The 
independence from S of magnetic reconnection is the conclusion of LV99 model. This 
model has been tested in Kowal et al. (2009, 2012) via a set of dedicated numerical 
simulations that confirmed the scaling predictions in LV99. This work exemplifies the 
advantages of the synergy of scaling arguments with numerical simulations as opposed 
to the "brute force numerical approach", which may not be productive while dealing 
with turbulence. 

One can argue that the process of reconnection diffusion was present in some of 
high resolution numerical simulations, although the researchers did not identify the 
process. However, without clear identification of the role of turbulence in fast recon- 
nection, one may not be sure when the results are due to physically motivated recon- 
nection diffusion and when they are the consequence of the bogus effects of numerical 
diffusion. For instance, Crutcher, Hakobian & Troland (2009) refer to the simulations 
in Luntilla et al. (2009) that produce, in agreement with observations, higher magne- 
tization of the cloud cores. If these cores are of the size of several grid units across, 
numerical effects rather than reconnection diffusion may be dominant and turbulence 
is suppressed at these scales. 

There are some very good news for numerical simulations, however. The LV99 
model predicts that the reconnection rates in turbulent fluids are independent of the lo- 
cal physics, but are determined by the turbulent motions. Therefore parasitic numerical 
effects that induce poorly controlled small scale diffusivity of magnetic field lines are 
not important on the scales of turbulent motions. In other words, the low resolution nu- 
merics may provide an adequate representation of high Lundquist number astrophysical 
turbulence as far as the reconnection diffusion is concerned. 

On the contrary, if the structures studied in numerical simulations (e.g. cores, 
filaments, shells) lose turbulence due to numerical diffusivity effects, we predict that 
the diffusion of magnetic flux in those simulated structures differs significantly from the 
diffusion in the actual interstellar structures where turbulence persists. As a result, naive 
convergence studies based on increasing the numerical resolution several times cannot 
notice the problem unless the resolution increased to the degree that the aforementioned 
structures become turbulent. 



9. Reconnection diffusion and simulations with ambipolar drift 

A number of studies attempted to understand the role of joint action of turbulence and 
ambipolar diffusion. For instance, Heitsch et al. (2004, henceforth HX04) performed 
2.5D simulations of turbulence with two-fluid code and examined the decorrelation of 
neutrals and magnetic field in the presence of turbulence. The study reported an en- 
hancement of diffusion rate compared to the ambipolar diffusion in a laminar fluid. 
HX04 correctly associated the enhancement with turbulence creating density gradients 
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that are being dissolved by ambipolar diffusion (see also Zweibel 2002). However, 
in 2.5D simulations of HX04 the numerical set-up artificially precluded reconnection 
from taking place as magnetic field was perpendicular to the plane of 2D mixing and 
therefore magnetic field lines were absolutely parallel to each other. This will not hap- 
pen in any realistic astrophysical situation where reconnection will be an essential part 
of the physical picture. Therefore we claim that a treatment of "turbulent ambipolar 
diffusion" without addressing the reconnection issue is of academic interest. 

Incidentally, the authors of HX04 reported an enhanced rate that is equal to the tur- 
bulent diffusion rate LVl, which is the result expected within the reconnection diffusion 
picture for the special set-up studied. The fact that ambipolar diffusion rate does not 
enter the result in HX04 suggests that ambipolar diffusion is irrelevant for the diffusion 
of matter in the presence of turbulence. This is another reason not to call the observed 
process "turbulent ambipolar diffusion'^. 

Therefore we believe that HX04 captured in their simulations a special degener- 
ate case of 2.5D turbulent diffusion where due to a special set up the reconnection is 
avoided and magnetic field lines do not intersect. We also note that, in the presence of 
turbulence, the independence of the gravitational collapse from the ambipolar diffusion 
rate was reported in numerical simulations by Balsara, Crutcher & Pouquet (2001). 



10. Summary 

1 . Magnetic reconnection and turbulence are inter-dependent processes. Magnetic 
reconnection gets fast in the presence of turbulence, while MHD turbulence requires 
magnetic reconnection to develop eddy-type structures. 

2. Turbulence in the presence of magnetic reconnection induces diffusion process 
that we term "reconnection diffusion" to stress the importance of magnetic reconnection 
for transporting magnetic fields and matter. 

3. Reconnection diffusion in most cases is the process that is faster than ambipolar 
diffusion. Therefore it is expected to dominate the removal of magnetic fields from the 
forming molecular clouds. 

4. Reconnection diffusion acts both in fully ionized and partially ionized gas. The 
effects of ambipolar diffusion are important only at the small scales and do not affect 
the important large scale dynamics of magnetic field diffusion. 

5. As diffusion of magnetic fields and matter do not depend on the microphysical 
processes turbulent volumes of numerical simulations represent the magnetic diffusion 
correctly in spite of tremendous differences in the Lundquist number of simulations and 
astrophysical reality. 

Acknowledgements.. The research is supported by the Center for Magnetic Self- 
Organization in Laboratory and Astrophysical Plasmas, the Vilas Associate Award and 
the NSF Grant AST-1212096. Hospitality of the International Institute of Physics (Na- 
tal, Brazil) is acknowledged. 



^ A similar process takes place in the case of molecular diffusivity in turbulent hydrodynamic flows. The 
result for the latter flows is well known: in the turbulent regime, molecular diffusivity is irrelevant for 
the turbulent transport. The process is not called therefore "turbulent diffusivity" without adding word 
"molecular". 
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